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Abstract 

A reduction theorem is proved for functionals of Gamma-correlated random fields 
with long-range dependence in d-dimensional space. In the particular case of a non¬ 
linear function of a chi-squared random field with Laguerre rank equal to one, we apply 
the Karhunen-Loeve expansion and the Fredholm determinant formula to obtain the 
characteristic function of its Rosenblatt-type limit distribution. When the Laguerre 
rank equals one and two, we obtain the multiple Wiener-Ito stochastic integral repre¬ 
sentation of the limit distribution. In both cases, an infinite series representation in 
terms of independent random variables is constructed for the limit random variables. 

Keywords: Hermite expansion, Laguerre expansion, multiple Wiener-Ito stochastic 
integrals, non-central limit results, reduction theorems, series expansions. 


1 Introduction 

This paper considers the family of Gamma-correlated random fields within the general 
class of Lancaster-Sarmanov random fields. Such a class includes non-Gaussian random 
fields with given marginal distributions and given covariance structure. The bivariate 
densities of these fields have diagonal expansions. Lancaster (1958) and Sarmanov (1963) 
idenpendently discovered these expansions in the context of Markov processes, namely, for 
dimension d = 1, and correlation function y(|x — y |) = exp (— c\x — y \), c > 0. This line 
of research was also continued by Wong and Thomas (1962), where Laguerre polynomials 
were used as well as Hermite and Jacoby polynomials, in Markovian settings. The exten¬ 
sion of these limit theorems, based on bilinear expansions, to the context of long-range 
dependent (LRD) processes was considered in Berman (1982,1984), and also for random 
fields by Leonenko (1999), Anh and Leonenko (1999) and Anh, Leonenko and Ruiz-Medina 
(2013), among others. That this class of random fields is not empty follows from the results 
by Joe (1997), who constructed the system of finite-dimensional distributions for a given 
bivariate distribution consistent with their marginal distributions, using the calculus of 
variations and the maximum entropy principle. Some properties of stationary sequences 
with bivariate densities having diagonal expansions, and their limit theorems were ob¬ 
tained by Gajek and Mielniczuk (1999) and Mielniczuk (2000). Specifically, in Gajek and 
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Mielniczuk (1999), long-range dependence sequences with exponential marginal 

distributions and its subordinated sequences are studied. In particular, processes of the 
form Zi = (X 2 + lf)/2, i = 1,2..., where {W}£7i and are independent copies 

of a zero-mean stationary Gaussian process with long-range dependence, are investigated. 
The asymptotic behaviour of a partial-sum process of the long-range-dependent sequence 
{G(Zj)}iT 1 , constructed by subordination from {Zi}^2 zl , is the same as that of the first 
nonvanishing term of its Laguerre expansion (see also Taqqu, 1975, 1979, in relation to 
central and noncentral limit theorems for long-range dependence processes in discrete 
time). In Mielniczuk (2000), different properties of bivariate densities (not necessarily 
associated with stochastic processes) are studied in the case where they admit a diagonal 
expansion, which is referred as Lancaster-Sarmanov expansion, including Mehler’s formula 
for bivariate Gaussian distributions, Myller-Lebedev or Hille-Hardy formula for bivari¬ 
ate Gamma distributions, among others (see, for example, Bateman and Erdelyi, 1953, 
Chapter 10). In particular, Mehler’s equality and Gebelein’s inequality are generalized. 
In addition, conditions are established for defining long-range dependence sequences satis¬ 
fying the reduction principle, by subordination to discrete time stationary processes. The 
present paper extends these results to the general setting of random fields with continuous 
d-dimensional parameter space, defined by a regular compact domain of M. d . In particu¬ 
lar, a reduction theorem is derived for Gamma-correlated random fields with long-range 
dependence. Some noncentral limit results are established for long-range dependence ran¬ 
dom fields constructed by subordination from chi-squared random fields, in the cases of 
function G having Laguerre rank equal to one and two. We will pursue in more details 
the chi-squared random field case where an explicit representation of the random field is 
available. 

The paper is organized as follows. In Section [2], we define the Lancaster-Sarmanov 
fields. In Section [31 we consider the case of Gamma and chi-squared random fields. In 
Section[H we prove the reduction principle for Gamma-correlated random fields. In Section 
[51 limit theorems are obtained for the case of functions of chi-squared random fields with 
Laguerre rank equal to one and two. We give a multiple Wiener-Ito stochastic integral 
representation of the limits. Infinite series representations of the limits obtained in Section 
[5] are obtained in Section [6j We establish infinite divisibility in Section [71 

2 The Lancaster-Sarmanov random fields 

We now introduce here the class of Lancaster-Sarmanov random fields with given one¬ 
dimensional marginal distributions and general covariance structure. Denote by £ 2 (^ 2 , J~, P) 
the Hilbert space of zero-mean second-order random variables defined on the complete 
probability space (D, J 7 , P). For a probability density function p on the interval (l, ?’), with 
—00 <1 < r < 00 , we consider the Hilbert space L 2 ((l,r), p{u)du) of equivalence classes 
of Lebesgue measurable functions h : (l,r) —> R satisfying 

h 2 (u ) p(u) du < 00 , p(u) > 0. 
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Let us also consider a complete orthonormal system {e k (it)}fcL 0 of functions in 
L 2 ((Z,r), p(u)du), that is, 


/ r 

e k {u ) e m (u) p{u)du = S k , m , (1) 

where S kiTn denotes the Kronecker delta function. We introduce the following condition: 


Condition AO Let {£(x), x E be a mean-square continuous zero-mean homogeneous 
isotropic random field with correlation function 

7(IMD = B(||x||) = Cov(£(0),£(x)), «E'. 

We assume that the densities 

P{u ) = -^P{£,(x) <u}, u e 

d' 2 

p(u,iu,||x-y||) = ^r^P{£(x) <u, £(y) <w}, (u,w) e (l,r) x (l,r), 

exist, and that the bilinear expansion 


p(u,w, ||x — y||) =p{u) p(w) 


1 + X]7 fc (ll x -y||) e k (u) e k (w) 


k =1 


( 2 ) 


holds, where 

OO 

X^ 2fc (ii x ii) < °°’ v ii x ii>°’ 

k=1 

and {e k (u)}^L 0 is, as before, a complete orthonornal system in the Hilbert space 
L 2 ((l } r),p(u)du). Assume also that eg(u) = 1. The symmetric kernel 


K(u, w, ||x - y ||) 


P(u,w, 11x - y||) 
p(u)p(w ) 


1 + “ y| \)e k (u)e k (w) 

fc=i 


(3) 


plays an important role. 

The series ([2]) converges in the mean-square sense if the integral 


I 2 = 



i Ji 


K 2 (u,w , ||x — y|| )p(u)p[w)dudw 

K 2 (u, w, ||x - y||)<iP{£(x) < u} dP{£( y) < w} < oo, 


where I 2 — 1 is known as the Pearson functional for the bivariate density p(u , w, ||x||) (see, 
for example, Lancaster, 1963). Then, the symmetric kernel K(u, w) belongs to the product 
space L 2 ((l , r ) x (l, r),p<8>p(u, w)dudw ) of square integrable functions on (Z, r) x (Z, r), with 
respect to the measure p <S> p(u, w)dudw. Thus, the kernel K defines an integral Hilbert- 
Schmidt operator on the space L 2 ((l,r),p(u)du). From the spectral theorem for compact 
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and self-adjoint operators (see, for example, Dautray and Lions, 1985, p.112), for each x, 
y £ R d , the kernel K admits the diagonal spectral expansion 


K(u,w) = p ( u '™' H x =^r fc (||x-y||)e fe (ti)e fc H, (4) 

p{u) p{w ) ^ 

where convergence holds in the space L 2 ((l, r) x ( l , r),p®p(u, w)dudw). Here, r k (||x — y||) 
is the sequence of eigenvalues, associated with the orthonormal system of eigenfunctions 
{e k (u)}^L 0 , which could also depend on x and y in a general setting. 

Thus, Condition AO postulates the expansion d4]) for the case where 

r*(x,y)=/(||x-y||), 


and e k {u) does not depend on x and y. Condition AO then implies 

/ r 

e k {u)p(u)du = 0, k> 1 

rr rr 

E[e n (£(x))e m (£(y))] = / / e n (u)e m (w)p(u, w, ||x - y\\)dudw 


e n {u)e m (w)p(u)p{w) 1 + ^ 7 fc (||x - y|| )e k {u)e k {w) dudw 


k =i 


= Sn t m 7 n (||x - yII), n, m > 1. 


(5) 


We will call the random fields satisfying Condition AO Lancaster-Sarmanov random 
fields, due to Lancaster (1958) and Sarmanov (1963). In the next section, we will refer 
to the special case of Gamma-correlated random fields, and, in particular, to the case of 
chi-squared random fields. We will also let (l,r) in (jT]) be (0, oo). 


3 Gamma-correlated random fields 

In this paper all random fields considered are assumed to be measurable and mean-square 
continuous. We refer to the class of random fields with Gamma marginal distribution 
and given correlation function. For details see Berman (1982,1984), Leonenko (1999), 
Anh, Leonenko and Ruiz-Medina (2013), among others. Following the ideas of Lancaster 
(1958) and Sarmanov (1963), we introduce a homogeneous and isotropic random field 
{£(x), x € R rf }, with given one-dimensional Gamma distributions, and given correlation 
structure 7 (||x — y||) = Corr (£(x), £(y)), x, y € R d . Let 

ppiu) = f^y u/3_1 ex P(-^)> u > °, P > ( 6 ) 

be a Gamma density, and let L 2 ((0, oo ),pp(u)du) be the Hilbert space of square integrable 
functions with respect to the measure pp(u)du, i.e., the space of functions F such that 

F 2 {u)pp(u)du < oo. (7) 
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An orthogonal basis of the Hilbert space L 2 ((0,oo),pp(u)du) can be constructed from 
generalized Laguerre polynomials Lp, k > 0, of index /3 (see Bateman and Erdelyi, 1953). 
Specifically, its elements are defined as follows: For k, m > 0, 


efc(«) = e k ] (u) = L[ p L \u) 


(yS-l)/ 


k\T(/3) 


1/2 


r°° 

, / e%\u) effi(u) pp{u)du = 5 k , m , (8) 

Jo 


T{P + k) 

where by Rodriguez formula for Laguerre polynomials 

= lP\u) = exp(u)—j: |exp(-M)tt /3+fc | . 

The first three polynomials are then given by 


(9) 


ef ?{u) = 1, e\ PI {u) = ^(p-u) 

eP\u) = {v?-2(J3 + l) u+ (0 + 1)0) [2 (/3 + 1) ^ . 


(/?), 


( 10 ) 


Applying Myller-Lebedev or Hille-Hardy formula (see Bateman and Erdelyi, 1953, 
Chapter 10) we obtain 


Pp(u,w, ||x — y||) = pp(u) pp(w) 


uw 


t(II x - y II) 


l + J]7 fc (ll x -y||) ef\u) ef\w) 

k =l 

C 8 -l )/2 


exp 


u + w 


x I/3-i 2 


1 -7(l|x-y||) 
1 


Vuw'y(\\x - yII) \ _ 

l-7(l|x-y||) / T(/3) (1 — 7(ll x — y||))' 


( 11 ) 


where y(||x — y||) is a continuous non-negative definite kernel on M d X M. d , depending on 
||x — y ||, and I s (z) is the modified Bessel function of the first kind of order g, with 

4(~) = t4? i 2 \ T\ [ (! - t 2 Y~ 1/2 exp (zt)dt, z> 0. 

Vvrr [g+ I) J_i 

Summarizing, one can define a homogeneous and isotropic gamma-correlated random 
field as a random field {£(x), x € K d } , such that its one dimensional densities 

Apk ( x)< u ] 

and two-dimensional densities 

8 2 

P{u,w, ||x - y||) = Q udw P [£(*-) ^ u ^(y) < w \ 
are dehned by © and m, respectively. In addition, the correlation function 7 satisfies 


X^7 2fc (ll z ll) < ll z ll>0- 

k =1 
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From equation (f7|) . F(u) can be expanded into the series 

OO 

F(u) = Y J c ^f\u) : C L q = 

q=0 

which converges in the Hilbert space L2((0,oo),pp(u)du). In particular, 

POO 

c o= F{u)ef- ] {u)pp(u) du = E[F(f(x))]. (13) 

Jo 

The Laguerre rank of the function F is defined as the smallest k > 1 such that 

Cf = 0, Ct^O. 

From equation m , for a homogeneous and isotropic Gamma-correlated random field 
{£(x), x G W 1 }, with correlation function 7 , the following identities hold: 

Ebf («(*))] = 0, E[e<®(£(x))4«(e(y))] = <W/(l|x - y||). (14) 

In order to introduce long-range dependence for Gamma-correlated random fields, we 
assume the following condition: 

Condition Al. The non-negative definite function 

7(l|z||) = m Z€M "’ 0 < 5 < d, (15) 

||z||° 

where C is a slowly varying function at infinity. 


F{u) ef\u)pp{u) du, q = 0 , 1 , 2 ,. 


( 12 ) 


3.1 The chi-squared random fields 

One can construct examples of random fields with marginal density ([ 6 ]) and bivariate 
probability density (1111) considering the class of chi-squared random fields. The chi-squared 
random fields are given by 

Xr( x ) = ^(U 2 ( x ) + --- + ^r 2 ( x )), xGM d , ( 16 ) 

where Yi(x),... ,Y r {pc) are independent copies of Gaussian random field {T(x), x G 
with covariance function F>(||x||) with F?(||0||) = 1. In this case 


7(ll x -y||) 


Cov(Xr( x ), Xr(y)) 
Var(y2(o)) 


H 2 (||x-y||), p = r/ 2 . 


(17) 


Note that by construction, the correlation function of chi-squared random fields is always 
non-negative. Moreover, 


EXr(x) 

and 


r 

2 ’ 


Vary 2 (x) = ^Var Y, 2 (x) = Cov(x 2 (0),Xr( x )) 

E[4 r/2) (Xr( x )) e^ /2) (Xr(y))] = <W B 2m ( ||x - y||), 



(18) 
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since as noted in dTJ) , 


coo 

/ e k /2 \ U ) e m 2 \ u ) Pr/2(u) du = S kjTn . 

Jo 

In the case of chi-squared random fields (1161) the analogous of Condition A1 setting 
in (1151) is the following Condition A2. 

Condition A2. The random field {Y(x), x € M d }, whose independent copies define the 
chi-squared random held (USD, is a measurable zero-mean Gaussian homogeneous and 
isotropic mean-square continuous random held on a probability space (Q,A,P), with 
EY 2 (x) = 1, for all x € M rf , and correlation function E[Y(x)Y(y)] = £?(||x — y||) of 
the form: 

B(\\z\\) = 1&T’ z€R "’ °<«< d / 2 - (!9) 

From Condition A2, the correlation function B of Y is continuous. It then follows that 
C(r) = 0(r a ), r —* 0. 


4 Reduction principle for Gamma-correlated random fields 

The following reduction principle is an analogous in spirit to the reduction principle of 
Taqqu (1975, 1979); for Gamma-correlated random helds, see also Berman (1982, 1984), 
Leonenko (1999), among others. 

From equation (fTP) . 


E 




4 /3} (^( x )) e m ) (?(y)) dxdy 


h,m °i( T )> 


( 20 ) 


where T> (T) denotes a homothetic transformation of a set DcR 11 with center at the point 
0 € T> and coefficient or scale factor T > 0. In addition, T> is assumed to be a regular 
compact domain, whose interior has positive Lebesgue measure, and with boundary hav¬ 
ing null Lebesgue measure. Dirichlet-regularity here is understood in the general setting 
established, for example, by Fuglede (2005, p. 253), as given in the following dehnition. 


Definition 1 For a connected bounded open domain T> with boundary dT> we say that 
xo € dT> is regular if and only if it has a Green kernel G v such that, for each x € T>, 

lim G®(x, y) = 0, Vy € V. (21) 

X—>X0 

The set T> is regular if every point of dT> is regular. 

Dirichlet regularity of domain T> ensures that the eigenvectors of the operator K a , 
introduced in equation (1371) below, vanish continuously in the boundary of domain T> (see, 
for example, Brelot, 1960, p. 137 and Theorem 32, in the context of potential theory, 
and, more recently, Chen et al., 2012, p.484, for 0 < a < 2, in the context of subordinate 
processes in domains). 
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In equation (1201) . under Condition Al, for 0 < 5 < d/k, 


a 2 k (T) = Var 


JV(T) 


7* (||x - y||) dxdy = [a d , k (V)\ 2 T 2d - kS C k {T){l + o(l)), (22) 


It>(t) Jv(t) 


as T —> oo, with 


a d ,kiP) = 



1 


| k5 


dxdy 


1 1/2 


, k > 1. 


/pip ||x - y| 

Note that, for the particular case of chi-squared random fields we have from (1181) 


(23) 


E 


!V(T) JV(T) 


e fcC\r( x )) e m(Xr(y)) dxdy 


— d k,m cr fc(^')) 


where, under Condition A2, for 0 < a < -M-, 


a 2 k (T) = Var 


lv(T) 


efc(Xr( x )) dx 


5 2fc (||x - y||)dxdy = [af k {V)fT za -^ a C Zk (T){ 1 + o(l)) 


2nn2d—2kot r*2k t 


IV{T) JV{T) 


as T —> oo, with 


= 



i 


\2kot 


dxdy 


1/2 


k > 1. 


LJx) Jd ll x - y 

The following theorem states the reduction principle. 


(24) 


(25) 


Theorem 1 Let |£(x), x E R d } be a Gamma-correlated random field. Assume that 
Condition Al holds, and that the function F E T 2 ((0, oo),pp(u)du) has generalized 
Laguerre rank equal to k, where pp(u) is defined by |ij). If the limiting distribution of the 
functional 


sk = 


a d , k {V)£ k / 2 (T)T d -(kS)/2 


F(£(x))dx-C 0 L T d |D| 


Id(t) 


(26) 


for 0 < 5 < d/k, exists as T —> oo, then it coincides with the limit distribution of the 
random variable 


ct 


dx - 


a d , k (V)C k V{T)T d -<~»- j vm 
The constants C k and C/f are defined in equations m and m, respectively. 


















Proof. The proof is based on the generalized Laguerre polynomial expansion of the 
function F. Specifically, under Condition Al, since 'y(||x||) < 1, and 7 ( 0 ) = 1, we have 


/c+Z / 


X 


< 7' 


k+l, 


l > 2 . 


Hence, from equation (l22l) . for T sufficiently large, 


a d +V)£ k / 2 (T)T d ~( kS )/ 2 


f F(£(x)) dx-C 0 L T d \V 
Jd(t) 



e f } (£( x )) dx 


< 


_a d:k (V)C k / 2 (T)T d ~( kS )/ 2 


OO 

V., (C > )2 L 


j=k+l 


V(T) JV(T) 


7 J (ll x -y||)^ x rfy < 


< 


1 2 


[a dik (l ))£t/ 2 ( T ) T i-(H )/2 j J V(T) J V(T) 


/ / 

JV(T) JV 1 


7 


fc+1. 


X - 


y||)cfocdy ^ {Cff = K R . 

j=k +1 


By Condition Al, for any e > 0, there exists Ao > 0, such that for||x — y|| > Aq , 
7(ll x — yII) < e- Let I<i = {(x,y) € V(T) : ||x - y|| < A 0 )}, and I< 2 = {(x,y) € V(T) : 
11x - y || > A 0 )}. Then, 


Jfc+i/ 


lv(T) JV(T) 


l x — y ll)^ xrf y = 


+ 


Ki 


K 2 


7 fc+ 1 (||x-y||)dxdy 


_ 0 ( 1 ) 1 0 ( 2 ) 

- kDrp kJrj-i 


(27) 


Using the bound 7 fc+1 (||x — y ||) < 1 on K\, and the bound 7 fc+1 (||x — y||) < e 7 fc (||x — y| 


on K 2 , we obtain, again, for T sufficiently large, 


S, 


(i) 


< 


J 7 fe+ 1 (ll x -y||)^y 


< Mi T d , 


for a suitable constant Mi > 0 , and 


S. 


( 2 ) 


< 


7 fc+ 1 (||x-y||)dxdy 


k 2 


< e 


[ [ 7 fc (||x-y||)dxdy 
J Jk 2 


< eM 2 T' 2d ~ kS £ k (T), 


for suitable M 2 > 0, and arbitrary e > 0. Thus, 

-1 2 r 


Kr = 


1 


_a d)k (V)£ k / 2 (T)T d -(kS )/2 
rpd 


7' 


k+l / 


< Ml V M 2 


It>(T) Jv(T) 

rji2d— k5 £k 


x-y||)dxdy J2 {Off 


j=k +1 


+ e- 


a l k {V)£ k (T)T 2d ~ kS al k (V)£ k (T)T 2d ~ kS 
which can be made arbitrary small together with e > 0 . 


(28) 


The following additional condition is assumed for the slowly varying function C in 
Theorem [2] below. 


9 




























Condition A3. Let C be the slowly varying function introduced in Condition A2. 
Assume that, for every m> 2 there exists a constant C > 0, such that 



f £(r||xi-x 2 ||) £(r||x 2 -x 3 ||) 

v A T )II X 1 - x 2 II^ >C(T)||x 2 - x 3 n* 5 


£(T||x m - xi||) 
£(T)||x m - X!II 5 


dxidx 2 • • • dx m < 


< C 


/ • • • (m) 

Jv 


dxidx 2 • • • dx r 


Jv ll x l - x 2||‘ 5 H x 2 - X 3 || 5 


X m - Xi 


Condition A3 is satisfied by slowly varying functions such that 


sup 

T,x i,X2GX> 


£(T|| Xl - x 2 ||) 

C(T) 


<C 0 , 


(29) 


for 0 < Co < 1. This condition holds, for example, for logarithmic type slowly varying func¬ 
tions £(||x||) = log(C + ||x||), C > 0, in the case where V C 13(0), with 

13(0) = {x e M d , ||x|| < 1}. 

Note that 


^(ll z ll) (1 4- || z ||/9)7 ’ ^ < P — 7 > 0; 

is a particular case of the family of covariance functions m studied here, satisfying 
Condition A3, with a = /Ly, and £(||z||) = ||z||^ 7 /(l + ||z||^) 7 . 

The next result involves chi-squared random fields. It provides the limit in distribution of 


1 

a^ 1 (P)£(T) T d ~ a 



Jr!*) 

e i 


(Xr( x )) dx > 


(30) 


and more generally, in view of Reduction Theorem [U of e)/ 2 being replaced by a function 
F with Laguerre rank 1. 


Theorem 2 Let (x 2 (x), x S M d } be the chi-squared random field introduced in ITU) , and 
consider the functional 


oXr 

FJrj-i 


a x d ^(V)C{T)Td- 


/ 

Jv(T) 


F(Xr( x ))dx - CjfT d \V 


(31) 


where ayfi (T>) is given in I25j) for k = 1. For 0 < a < d/2, under Conditions A2 and 

,2 

A3, its limit, in distribution sense, S&H , in the case of F having Laguerre rank k = 1, has 
characteristic function of the form 


<t>{z) = Eexp{izS%£} = exp (?- ^ ( 2i«/y/2r) g \ ^ z € M, 

\ 2 z ' m j 

where c m , m > 2, are defined as follows: 

c m = f ■ ■ ■ f t, - 7— 7 - 7~ ■ ■ ■ 7 -— clxi ... dx m 

Jv M Jv II X 1 - X 2 || a ||x 2 - X 3 || a ||x m - Xi||“ 


(32) 


(33) 
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Remark 1 Note that, from Theorem [7J applied to the particular case of chi-squared ran¬ 
dom fields with k = 1, 


oXr 

&OO 


= lim 

T 


Cf 


oo n X r 
a d, 1 


(' D)C(T)T d ~ a Jv(t) 


ef /2) (Xr( x )) dx - 


(34) 


_ yZ 

Proof. From Remark [T] (see equation m and Theorem [T]) , the limit distribution of SjT 
as T -> oo, if it exists, can be obtained as the limit in distribution given in (1341) . since F 
has Laguerre rank k equal to one. 

The first Laguerre polynomial of the chi-square random field (x), x £ M. d } is the 
sum of r independent copies of the second Hermite polynomial of the original Gaussian 
random field (T(x), x £ M rf } involved, that is, for x £ M d , we have by (1101) . 


e! r/2) (Xr( x )) = 


(2 E y /( x )l ^EEE) 1} V^^ H2(Yj 

\ 3 =1 ) 3 =1 3 = 1 


2 / r 
r 


( x ))- 

(35) 

From equation (1351) . one can prove, in a similar way to Theorem 3.2 by Leonenko, Ruiz- 
Medina and Taqqu (2014), that the limit characteristic function admits the expansion 
Specifically, 


4>t 0 ) = E 


iz 


exp 


n-P 

3= 1 


T d-a£( T ) a Xr( P ) 7 x,( T ) 

1 ^ 1 / -2 iz 


dx 


E 


= 2 m \\/ 2 ra*\(V)T d ~ a £(T) 


Tr UVP MT) 


= exp 


-E 1 f- 

V 


—2 iz 


y/2Ta^(V)T d ~ a C(T) 

Note that under Condition A2, since EY 2 (x) = 1, 


Tr ( R^ v(t) 


(36) 


dx = 


IV(T) 


/ E [y 2 (x)] dx = E 
J V(T) 

00 

3 =1 


!v(T) 


Y 2 (x)dx 


E \j,T{Ry,v(T))En 

3 =1 


In the study of the convergence of the series (1361) . to apply Dominated Convergence 
Theorem, we use Theorem 3.1 by Leonenko, Ruiz-Medina and Taqqu (2014), where it is 
proved that, for 0 < a < d/2, the squared K? a of the operator 

£a(/)(x) = [ |T- ] —^f(y)dy, V/ £ Supp (JC a ), 0 < a < d, (37) 

Jv ll x — yii 

is in the trace class. In particular, its trace is given by 

Tl ' (E) = [ I .I 1 < ]2a dxd y = \°u/\ ( V )} 2 < oo- (38) 

Jv Jv ii x — y 11 
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From the Definition of the Fredholm determinant of a trace operator (see, for example, 
Simon, 2005, Chapter 5, pp.47-48, equation (5.12)) the Fredholm determinant of K? a is 
given by 


D /Ci(u) = det(I - uK? a ) = exp ( - jT 


TrJC 2k 


k 


= 


(39) 


for u) € C, and |w|||/Cq||i < 1, with ||/C^||i denoting the trace norm of operator K? a . In 
particular, for u = 2iz, and for \z\ < ^ , 


[D Kl {2\z))- 1 / 2 = exp (j ^(2iz) fc j < oo. 

In addition, under Condition A3, there exists a positive constant C such that 


(40) 


^2 Tr ( R Y,V(T) 


V JV 


£(T||xi-x 2 ||)£(r||x 2 -xi||) 1 

2a 


c 



C(T) 

1 


V JV ll x l — x 2 


C(T) ||x 1 -x 2 | 

j 2 ^dxi(fx 2 = CTr (/Cf) < oo, 


dxidx 2 


(41) 


Tr ( iZ™ m 
d m V Y ’ V ( T ) 

i r 


£(T||xi - x 2 1|) £(T||x 2 - x 3 1|) C(T ||x m - xi| 


< 


[C(T)} m J V (m) J V || x l x 2 || a 

1 


| x 2 - x 3 | 


| x m - x l| 


rixi ... dx n 


cf ' 

( m ) Jv Il x i — 


x 2 || a ||x 2 — x 3 || q 


| x m - x l| 


dxi . . . dXr, 


= CTr(K%) < oo, m > 2, (42) 

since ||/C™||i < jq/C* || 1; K > 0, for m > 2. Here, d T = a£f 1 (P)T d - a £(T). 


From equations (|36l) and dT0l) - (H2l) . for 0 < z < y/r/2 A Jr/ (2||/C^, ||f), we obtain 


I </>t( z )\ < 


rC ^ i 


exp 


E 

=2 

oo 


2 ' m 

m=2 


rC 


exp - 


< 


2 r) TV(/C-) 

E— ' 

' 2m 

771=1 

2m + 1 

OO 

^-(-2U/v^)"' T r(C5”) 
2r)' Tr (/Cf, m ) 


_ x 2m _ 

2r) Tr(JC 2 a m ) 


, -X 2m+l 

2r) Tr (JC 2 a m+1 ) 


m= 1 


rC 


exp 


_m= 1 


+ V — 

z — J 771 
m=l 


“ Va/21 


-rC 


< oo, (43) 
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where the last identity in (14311 is obtained from the definition of the Fredholm determinant 
oiKl as given in equation (1391) . 

We can thus apply the Dominated Convergence Theorem to obtain lim t^>-oo' 1 Pt( z ) = 
ip(z), for 0 < z < \]r/2 A yjr/ (2||/C„||f). An analytic continuation argument (see Lukacs, 
1970, Th. 7.1.1) guarantees that if defines the unique limit characteristic function for all 
real values of 2 . 


5 Limit theorems for Laguerre rank equal to one and two 
and Wiener-Ito stochastic integral representations 


Consider the chi-squared field defined in (1161) . The multiple Wiener-Ito stochastic integral 
representation of the limit in distribution of the functional (131 D . and of the functional 


S r 2,T 


1 

af 2 (V)C 2 (T)T d ~ 2a 



,(r/2) 


(Xr( x )) d* 


(44) 


is derived in Theorems|3]and|l]below, respectively. Here, afff 2 is defined as in (|25l) for k = 2, 

and C(T) is the slowly varying function introduced in (fTUl) . The basis function e 2 2 {u) is 

defined in Relation ([21) • In the Section [6J from the multiple Wiener-Ito stochastic integral 

representations derived in this section, we obtain an infinite series representation, in terms 

,2 

of independent random variables, for SSS, with characteristic function given in Theorem 
El and for the random variable obtained as the limit in distribution of the functional (1441) 
in Theorem E] below. 

The slowly varying function C in ()19[) is assumed to belong to the class CC which is 
now introduced (see Definition 9 by Leonenko and Olenko, 2013). 


Definition 2 An infinitely differentiable function £(•) belongs to the class CC if 

1 . for any 5 > 0, there exists Ao(5) > 0 such that \~ s C(\) is decreasing and X s C(X) is 
increasing if A > Ao(<5); 

2. Cj G SC, for all j > 0, where Cq(X) := C, £j + i(A) := A£' (A), with SC being the class 
of functions that are slowly varying at infinity and bounded on each finite interval. 


The following lemma will be applied in the proofs of Theorem El and [4] below (see 
Theorem 11 by Leonenko and Olenko, 2013). 

Lemma 1 Let a G (0, d), S G C'°°(s ri _i(l)), and C G CC. Let (A(x), x G K d } be a mean- 
square continuous homogeneous random field with zero mean. Let the field X has spectral 
density /( u), u G M d , which is infinitely differentiable for all u / 0 .If the covariance 
function -B(x), x G M. d , of the field X has the following behavior 

( a ) ll x ir^( x ) ~ S (pjf) £(ll x ll), x — >°o, 

the spectral density satisfies the condition 
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(b) ||u|| d - a /(u) ~ S a , d (^) C (^) , ||u|| —► 0. 

In Propositions [Tj [2] and Theorems EHH below, the following Fourier transforms and con¬ 
volution formulae will be applied in S(R d ), the space of infinitely differentiable functions 
on M d , whose derivatives remain bounded when multiplied by polynomials, i.e., whose 
derivatives are rapidly decreasing (see Lemma 1 of Stein, 1970, p.117). 


Lemma 2 (%) The Fourier transform of the function ||z|| d+/d is z/(/3)||z|| in the 

sense that 

f \\z\\~ d+>s ip(z)dz = f u(P)\\z\\~^ F(ip)(z)dz, \/ip € <S(M d ), (45) 

J R d Jm. d 


where 


and 


v(P) = i 0 < p < d, 


d-p 

2 


(46) 


F(ip)(z) = / exp(—i(x,z ))'0(x)dx 

Jm d 

denotes the Fourier transform of ip. 

(ii) The identity F ((—A )~^^ 2 (f)) (z) = ||z|| _ ^J r (/)(z) holds in the sense that 

[ (-A)' /3/2 (/)(x) 5 (x)dx= —^ f J r (/)(x)||x||~' 3 7 r ( 5 )(x)dx, Vf,geS(R d ), 
jR d (2ir) a J R d 

(47) 

for 0 < /3 < d. 

(in) The following convolution formula is obtained by iteration of [J7\ ) 


L 


i 


R d z/(4/3) 


\~ d+ip f{z)dz = / ||z||- 4 ^(/)(z)dz 


Im d W)] A Im 


|z - xi||- d+/3 ||xi - x 2 ||-^||x 2 - y11 ~ d +P 11 y 11 ~ d +P ^xidx 2 dy 


x/(z)dz, V/ £ 5(1“), 0 < /3 < d/4. 


(48) 


The proof of this lemma can be seen in Stein (1970, p.117), and Leonenko, Ruiz-Medina 
and Taqqu (2014). 


Proposition 1 For 0 < a < d/2, the following identities hold: 


[ \K(X 1 + X 2 ,V)\ 2 —- 

JR 2d (||Ai 


dX\dX2 


\d—a 


a ii(p) v ( a ) 


-i 2 


|P| 


Ha)] 2 Tr(JCi) 

\V \ 2 


< oo, 
(49) 
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where a/fi (V) is defined in f%5l) . v(a) is introduced in equation and K is the char¬ 

acteristic function of the uniform distribution over set T>, given by 

K (X,V) = [ e-'^p v (x) dx = -L [ e-’^dx = ^, (50) 

with associated probability density function px> (x) = 1/ \V\ if x eP, and 0 otherwise. 

Remark 2 Note that for D = £3(0) = {x G ||x|| < 1}, the function i?(A) in \5U\) is of 
the form: 

f exp (i (x, A)) dx = (27r) d/2 ^ ||]^^ , d > 2, 

Jb( o) WM ' 

where J v {z) is the Bessel function of the first kind and order u > —1/2. For a rectangle, 
F> = [I = {ai < Xi < bi, i = 1,... ,d} , 0 G n< 

d 

m = n (exp (iA jbj) — exp (iA jafi) /iA j, d > 1. 

3 = 1 

(see, for example, Leonenko and Olenko, 2014)- 

Theorem 3 Assume that Conditions A2-A3 hold, 0 < a < d/2, and that C G CC. 

2 

Consider Sf r be the functional m, given in terms of the integral of functional F of 

the chi-squared random field with Laguerre rank equal to one. As T —>■ oo, the limiting 
2 2 

distribution of SfT , with characteristic function H32\) , admits the following double 
Wiener-Ito stochastic integral representation: 


qXr _ 
&oo — 



Zj (dX- L ) Zj (dX 2 ) 

V 1’ 2 ) d—a d—ac 

||Ai || 2 11A 2 11 2 


(51) 


where Zj, j = l,...,r, are independent Gaussian white noise measures, v is defined in 
m and the notation f^ 2 d means that one does not integrate on the hyperdiagonals X± = 
iA 2 . Here, 

H (A,, A 2 ) = K (Ai + A 2 , T>), (52) 

where K (X,V) is defined in LhOll . 

The proofs of Proposition |T| and Theorem [3] can be derived as in Theorem 4.1 by 
Leonenko, Ruiz-Medina and Taqqu (2014), since from Theorems [T] and El Theorem O 
holds for the functional S\_t given by (see (1M1L 


Si,T = 


a%\(p)C(X)T d ~ a Jv(T) 


e!' / 2 ) (Xr( x )) dx 


a^{V)C{T)T d ~ 


V2f 


E 


= 1 JV(T) 


H 2 (Yj(x))dx 


(53) 


We now turn to the case k = 2. 
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Proposition 2 Let T> be a regular compact set and let K (A, T>) be defined in H50\) . 
For 0 < a < d/4, the following identities hold: 



+ A 2 + A 3 + A 4 , P )| 2 


rui ( ^> 


nti diAi 


\d—a 




\v\ 


< 00 , 


(54) 


where a/f 2 (V) is defined as in equation f£5j) for k = 2 , and n(a) is introduced in equation 

m- ’ 

Proof. The proof follows from the application of Theorem 3.1 in Leonenko, Ruiz-Medina 
and Taqqu (2014), where the asymptotic spectral properties of operator KL a in equation 
(GZD, on a Dirichlet regular compact domain T >, are established. Let us now consider the 
following norm on S(R d ), 


II/IIU’--/-= (<" A 

= f (-Af"-« ! (/)( x )7w* = f ,. 1 . - / i,—L^/(y)7w<iy<ix 

Jud J R d v(d - 4a) J R d ||x - y || 4 

= TTFd! l^(/)(A)| 2 ||A||-( d - 4 “)dA, V/eS(M d ), 0 < a < d/4. (55) 

(2vr) a J R d 

The space LL^a-d = 5(M d / ^(-A ) 2 “- d / 2 j s th e jjilbert space of the functions of 5(M d ) with 
the inner product 

(/,d)(_A ) 2 «-<*/ 2 = [ 7 , 1 7 T [ |T-(y)fl'( x )^y^x, V/, 5 € 5 (M d ), (56) 

^ ’ JRd u{d - 4a) J R d ||x - y || 4 

and the associated norm (1551) . Here, S(R d ) II II( a) 2 “ d / 2 denote the closure of S(M. d ) with 
the norm (l55l) . Note that Equations (1551) and (l56l) can be extended to the space Ti^a-d by 
continuity of the norm. In particular, 


11x51 


1 


^4a — d 


1 ,, [«& m 2 

-r^dydx = - 


lx> v{d — 4a) Jx> ||x — y|| 4a u(d — 4a) 

As noted before, from Theorem 3.1 by Leonenko, Ruiz-Medina and Taqqu (2014), 

1 


(57) 


Tr (ICi) = 


Thus, for a = 2/3, 


Therefore, 



'v Jv ll x ~y 
1 


1 2a 


dydx < 00 , 0 < a < d/2. 


(58) 


vJv ll x -y|| 4/3 


dydx < 00 , 0 < /3 < d/4. 


[a X df 2 m 2 = 


1 


14a 


dydx = 1 fid — 4a)\\li>\\‘ii ia _ d < 00 , 0 < a < d/4. 


IvJv n x -y|| 

Equivalently, 1-p belongs to the Hilbert space Ti^a-d, for 0 < a < d/4. 
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Applying the convolution formula (14811 in Lemma [2] we then obtain 


feW _ lh l|2 \V\ 2 f ljr , _. |2|l 

j/(d-4a) 11 (2vr) rf 7 R J ^ l|cJl 

|T>| 2 i/(4a) 


|- d+4a du>i 


(2vr) d Ka)] 4 V 




, , , , ||— d+a\\, . , , ||— d+a 

UJl — U>2 \\0J2 ^3 || 


vll, . , , ||—'d+a 11 . , 11 — d+a 

X | Cc^3 — U 74 U>4 




Hence, 


|P| 2 t/(4a) f 
(2n) d [i'(a)] 4 J W 4d 


K>>' 2 - iS 


i=2 

2 


du;i 


WE A *’ P 


v i=l 


ru 

uU\\\\\ d - a 


K (E A ' /D 


V Z=1 


n 


n»=i 

4 "AJI^-" 1 


i= 1 IK'U 


since 


u( laMdd la ) _ ^ Equation (1M1) then holds. 


— Haql—cL) 


Note that by continuity of the norm in H^ a -d 

lx>*lz>(x)= f lx>(y)lx>(x + y)dy = f 1d(x + y)dy G L 2 (V) 

J R d Jv 

since 

/ / lx>(y)lx>( x + y)dy dx < |H i j (X)) (0)| 3 , 

JR d JR 11 

where |Z5 _r(x>) (0)I denotes the Lebesgue measure of the ball of center 0 and radius R(T>), 
with R{T >) being equal to two times the diameter of the regular compact set T> containing 
the point 0. Hence, F{lx> * lr>)(A) = \'D\ 2 \K (A, V)\ 2 belongs to the space of Fourier 
transforms of functions in R^ a _ ( j / (see also Remark 3.1 by Leonenko, Ruiz-Medina and 
Taqqu, 2014). 

■ 

Theorem [2] provided the limit of (13011 involving e ^ 2 . The next theorem provides the 
limit of flSD, involving e?/ 2 . Note that el/ 2 is defined in (fTOT) . but also satisfies (l6TT) below. 


Theorem 4 Assume that Conditions A2-A3 hold, and that C € CC. Then, for 
0 < a < d/4, the functional 5 2 ,t defined in 0 converges in distribution to the random 
variable admitting the following multiple Wiener-Ito stochastic integral representation: 


Soo = 


|X>| 


d 4[i/(a)] 5 


r( 2 + 1 


I- 1/2 


E f f k(±kv 

jfatj J R J R V 

k=l’ R Vi=l / n*=l 


Zj(d\i)Zj{d\2)Zk(d\^)Zk(d\i) 

nt 1 iiA l ||(^>/2 


nE z k {d\j) 

l l Ay || (^ — a )/2 


(59) 
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where the random measures Zj(-), j = 1,2,3,4, are independent Wiener measures. K(-,T>) 
is the characteristic function of the uniform distribution over the set V. The stochastic 
integrals f R2 d appearing in the first sum of L59\) are defined as mean square integrals, in 
which integration is excluded over hyperdiagonals X\ = ±A 2 , and A 3 = ±A 4 , related to 
each component Zj and Z k (see Fox and Taqqu, 1985). In the second sum, f^ 4d means 
that one can not integrate on the hyperdiagonals A, = ±Aj, i 7 ^ j, i,j = 1,2,3,4. 

Proof. 

The restriction to V(T) of the independent copies Yj, j = 1,..., r, of Gaussian random 
field Y, i.e., {Y)(x), x € T>(T), j = 1,... , r}, satisfying Conditions A2—A3, admit the 
following stochastic integral representation: 

y l( x ) = 7TTT / exp(i(x,A))A(A,P(T))/ 0 1 / 2 (A)^(dA), x g V(T), j = l,...,r. 

JR d 

( 6 °) 

It is well-known (see, for example, Anh and Leonenko, 1999) that 

4 r/ 2 ) (Xr( x )) = ^( r (^ + 1 )) ' 

r r 

x Y H 2 (Y k (x))H 2 (Xj(x)) - J2 ^4(ifc(x)) , 

k,j= 1 , k^j k= 1 

(61) 

where, as before, ,\'^(x) is the chi-squared random field introduced in (usd, and e' 2 ^ 
denotes the second Laguerre polynomial with index r/2 (see Bateman and Erdelyi, 1953, 
Chapter 10). Here, H 2 (u ) = u 2 — 1 is the second Chebyshev-Hermite polynomial, and 
H±(u) = u 4 — 6 u 2 + 3 is the fourth Chebyshev-Hermite polynomial. 

From equation (1611) . the functional (1441) admits the following representation: 


S 2 .T 


1 

af 2 (V)C 2 {T)T d - 2 “ 



4 / 2 (Xr( x )) dx 



- 1/2 1 
dx 


Y [ H 2 (Y k {x))H 2 (Y j {x))dx 

k,j= 1 , 


-Yf H 4 (F fe (x))dx 


(62) 


Using Ito’s formula (see, for example, Dobrushin and Major, 1979; Major, 1981), we 
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obtain from equation (16211 


o 1 / (r \\-V2 1 

S ^ T = 4 ( r \2^ )) T t 


V / 

fcj=l, fc#r® (T > JR2<i JR2d 


eX P \ X ! X! A * 


2—1 


si 


r*// 


/ exp < x 

“i ■'©co ■ /R4d V \ i=i / / r=i 


n yMran n 

2 — 1 2 — 1 2=3 

,e 4 nvww 


iLfr + or 172 w 

4 V \2 // [^(a)] 2 ^ 


E 



k \L i x " v ) 

x ([■'(a)] 2 n v/mmto ) n n z ‘< dA <) 


2=1 


/ 2=1 


2=3 


e £ 4d K (e a - p ) . 

/c —1 \ 2 — 1 / \ 2—1 / _ 


(63) 


Hence, applying Minkowski inequality, 


E 



< 


'1 

Ky>)] 

-1/2 |p| 

4 

[k«)] 2 J 

'1 

K£+0] 

-1/2 \V\ 

4 

K«)] 2 . 


E [Y 1T - Y, + y 2 - Yzr]' 


£ [Hit - Hi] 2 ) V ' + (e [Y 2 - Y 2T ] 2X ^ 


(64) 
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where 


E[Y xt - Yi] 2 = Y, E 

k,j;k^j 


f f We*.®) ^iiv®® 

J R 2dJ R 2d J [ d T fJl 


Zj{d\i)Zj{d\2)Zk{d\^)Zk{d\/C) 


- £ 

r 

E[Y 2 -Y 2T } 2 = Y. e 

k= 1 


nti iiA,n<‘ i -">/ 2 


K ’(£ Ai ’ I> 


V, 2=1 

4 


<2t(Ai, a 2 , a 3 , a 4 ) ■ 


d\d 


i =1 


|A ? ;|K- C 


[ k ( y ,\, v ) M^_n/wraM- 
Vs / L * r=i 

xZfc{dAi)Zfc(dA2)2ft(dA3)Z^(dA4}|‘ 


ns imi^ 2 


£ 

fc=i • 




V 2=1 


Qt(Ai, A 2 , A 3 , A 4 ) J ] 


dK 


i =1 


I A,- || rf —« : 


(65) 


with 


2 4 


«r=(n iia,ii i^ifLnvVwram -1 ]. 


dT 

The convergence to zero of Qt, as T —» oo, can be proved as in Theorem 4 .1 (ii) in 
Leonenko, Ruiz-Medina and Taqqu (2014) (see also Leonenko and Olenko, 2013). Hence, 
from equations (1M1) and (f65|) . as T —>• oo, 


E 


1 


S 2 ,T — T 


4 L \2 


r - + 1 


-i /2 mi 


K«)F 




4 . j./A | (V/A.)A, !’ 4 l./A I 

n‘„ 


>V 2 = 1 

' 4 


£ /.T £ A -® 


fc=l ’ 


v 2=1 


n,Lzt(^,) 

n,Li iiA.ii <' i -“)/ 2 


converges to zero, which implies the convergence in probability and hence, in distribution 
sense, as we wanted to prove. ■ 


6 Series representation in terms of independent random 
variables 

V 2 

We provide here series representations for the limit random variable SsS obtained when 
the Laguerre rank equals one, and for the random variable Soo obtained when the Laguerre 
rank equals two. 

Theorem 5 Assume that the conditions of Propositions QUH and Theorems & and\4\ hold. 
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(i) For the case of Laguerre rank equal to one, the limit random variable S$S in Theorem 
0 admits the following series representation: 


QXr _ 

^OO — 


-J I I 

-T7T m E E -!) = E E 


e 2 1 n -l), ( 66 ) 


where n(a) is given in { £ j n , n > 1 , j = 1 ,... ,r} are independent and identi¬ 

cally distributed standard Gaussian random variables, and 


A n (5&) = - 


y/2rv(a) 




with p n (TL), n> 1, feeing a decreasing sequence of non-negative real numbers, which 
are the eigenvalues of the self-adjoint Hilbert-Schmidt operator 


H(h){ Ai) = / i/i(A 1 -A 2 )/ l (A 2 )G Q (dA 2 ) : L 2 ^ 


L G a 


being 


G a (dx) = 


1 


X 


Id—a 


dx. 


(67) 


( 68 ) 


7/ere, f/ie symmetric kernel H\ (Ai — A 2 ) = H( Ai, A 2 ) = Ji , with K 

being as before the characteristic function of the uniform distribution over the setT>. 

(ii) For the case of Laguerre rank equal to two, the limit random variable Soo in Theorem 
[ 7 ] admits the following series representation: 


S n 


oo oo oo 


4Ka)P V (i + 1 )] 


- 1/2 


v\ 


EEE^ nWlpnl qn 


n= 1p=lg=l 


E “ 1 )( £ fc,«,n “ !) - E(^.» “ 1 )( £ fc,«,n ~ ^ 

k,j:k^j 


k =1 


, (69) 


where {£j, p , n , j = 1, ■ ■ ■, r, p > 1, n > 1} are independent standard Gaussian ran¬ 
dom variables, in particular, E[e j:P:n £ k , q ,m] = 5 n , m hp, q 5j,k, for every j,k = 1 ,..., r, 
and n,m,q,p > 1 . 7/ere, /i n (?7), n > 1, are f/ie eigenvalues, arranged in decreas¬ 
ing order of their modulus magnitude, associated with the eigenvectors ip n , n > 1 , 
of the integral operator H : L Ga ^ Ga (M 2d ) —>• L Ga ^ Ga (M 2d ) given by, for all 


h G L G a ®Ga 


p2 d 


)> 


77(d)(A 


.iA 2 ) = / K (EA 4 ,P ] d(A3,A 4 )G a (dA 3 )G Q (dA4). 


(70) 


Additionally, for each n > 1, 7 jn , j > 1, are ide eigenvalues, arranged in decreasing 
order of their modulus magnitude, associated with the integral operator on L G (M rf ) 
defined by kernel 
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Proof. The proof of (i) can be derived as in Corollary 4.1 (see Appendix A) in Leonenko, 
Ruiz-Medina and Taqqu (2014). 

(ii) From Proposition [2l the operator defined in (1701) is a Hilbert-Schmidt operator. Equiv¬ 
alently, the kernel 

H(X 1 ,X 2 ,X 3 ,\ 4 ) = k(^2 X ^ V 

V 2=1 

belongs to the space L 2 Ga ®G a ® Ga ®G a ( R4d ) ■ Thus > ^ € 5 ( L G a ®G a ( rM )) > where as usual 
S(H ) denotes the Hilbert space of Hilbert Schmidt operators on the Hilbert space H. 
Hence, it admits a kernel spectral representation in terms of a sequence of eigenfunctions 
{tp n , ^ > 1} C L 2 Ga ^ Ga (M 2d ) , and a sequence of associated eigenvalues {p n (T~L), n > 1}. 
That is, 

OO 

H(X i, A 2 ,A 3 , A 4 ) = fi n (W.)tp n (\i, A 2 )y> n (A 3 , A4). (71) 

n =1 

In particular, since, for every n > 1, p n € L Ga ^ Ga (M m ) , then, 



a 2 )| 2 


Ai 


d\\d\2 



< OO, 


which means that <^ n (Ai, A 2 ) defines an integral Hilbert-Schmidt operator T on Lg, (M d ), 
given by 

T(/)(A 1 )=/ <p n (Xi,X 2 )f(X 2 )G a (dX 2 ), Vf e L 2 Ga (R d ). 

Therefore, it admits a spectral kernel representation in L 2 G (M rf ), in terms of a sequence 
of eigenvalues {7 pnl p > 1 }, and an orthonormal system of eigenfunctions {4> pn , p > 1 } of 
Lg a (M d ), of the form 


7 ? n(Al)A 2 ) ^ ] '7pn < ^pn(Al)0p n (A 2 ), (72) 

i G a ®G a ( RM ) p=i 

for each n > 1, where convergence holds in the norm of the space (M 2rf ). Replacing 

in equation (17T1) the functions {<^ n }^ =1 by their respective series representations as given 
in equation (1721) . we obtain 


OO OO OO 

3,a 4 ) 2 = X£I> nftty'l pn'Y qn$pn(X 1 ) 4 > pn{X2^) 4 > qn{X3^) 4^ qn .(A 4 ), 

^ n=lp=lq=l 

(73) 


where convergence holds in the norm of the space L 2 &iG (M 4<i ) := L Ga ^ Ga( % Ga( %> Ga 
since, from equations (1711) (1721) . considering Minkowski inequality, we have 
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(74) 


where we have applied convergence in Lg ag)Ga (M 2d ) of the series 7 qn&qni') ® rf’qni') 

to the function tp n (-,-), for each n > 1 , which, in particular, implies that such a series 
differs from <p n (-,-) in a set of null G a < 8 > G a -measure. 

From Theorem |4l 
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Replacing #(Ai, A 2 , A 3 , A 4 ) = iF A*, £>) by its series representation (1751) in the 

above equation, one obtains 
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Applying Ito’s formula (see, for example, Dobrushin and Major, 1979; Major, 1981), 
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(75) 


as we wanted to prove. 

In addition, the orthonormality of the systems of eigenfunctions n > 1}, in the 
Hilbert space 7v Ga 0 Ga (M 2d ), means that 


{.^Pni ^Pk) l% „ (R24) — / ^(Ai, A 2 )</7 fe (Ai, A 2 )G a ((7Ai)G a ((iA2) — <5 n fe, (76) 

G a ®G a V. ' J R 2d 

where, as before 4 n ./- denotes the Kronecker delta function. Replacing (p n and </? fc in (1761) 
by its series representation (1721) in L Ga (R d ), we obtain 
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which implies that 
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and we know that 


n = k, 


(79) 


(fipm&qk) L^(Kd) ^P,9> 

from the orthonormality of the system of eigenfunctions providing the diagonal spectral 
representation (|72l) of <p n , for each n > 1. Hence, in addition, from (1771) and ([79]) . we have 

OO 

X^7pn = 1 > Vn>l. 

p=i 

Thus, from equations (1771) (1771) . and from the independence of the Gaussian copies YJ(-), 
* = 1, ■ ■ •, r, of random held Y (•), we obtain T[£j, ? , n £jfc,p,m] = S n , m S Ptq 6 k}j , for every j, k = 
1 ,... ,r, and n, m,q,p > 1 . 


Corollary 1 Under the conditions of Theorem 0 for Laguerre rank equal to two, 
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(80) 


where l 1 and 1 respectively are 1 X r and r x 1 vectors with entries equal to one. For 
p,n> 1, £p. n denotes a r x 1 random vector with entries e 2 n — 1, / = 1 ,..., r, ( 8 ) denotes 
the tensorial product of vectors, and Trace(A) the trace of a matrix A. Thus, S 0 0 admits 
an infinite series representation in terms of the sequence of independent random variables 
{r] n , n > 1 } given in [8(h) . 

The {rj n , n > 1} are independent because, for each n > 1, is a function of random 
variables {(e 2 n — 1), i = 1,..., r, p > 1} and, as follows from their definition in equation 
(USD, for n / k, with n,k > 1, {(e 2 p>n - 1), i = 1,..., r, p > 1} and {(e 2 g fe - 1), i = 
1 ,...,r, q > 1} are mutually independent, since the function sequences {4 > pn } P >i and 
{4> q k}q>i ar e orthogonal in the space Lg, Q (M d ) (see equation (jTHD ). 


7 Infinite divisibility 

y 2 

When the Laguerre rank equals one, Soo is infinitely divisible. Theorem [5](i) allows the 
derivation of its Levy-Khintchine representation. It is given by 

Theorem 6 Under the conditions of Proposition [7] and Theorem 0, 


<f>(9) = E 


exp (i 9S& 


= exp 


(exp(iud) - 1 - mO) p a/d (du) , 


(81) 
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where p a /d is supported on (0, oo) having density 


= £ £ 


exp 
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, u > 0. 


(82) 


Furthermore, q a /d has the following asymptotics as u —> 0 + and u —> oo, 
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where 


c(d, a) = it' 


= W 2 
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r(§) [ r (j)] <d_ “ )/ ‘ i 

Proof. Let us first consider a truncated version of the random series representation (1661) 

r M 


1=1 k =1 


with S££ —» <S2b r , as M tends to infinity. From the Levy-Khintchine representation of the 
d 

chi-square distribution (see, for instance, Applebaum, 2004, Example 1.3.22), 
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(84) 


To apply the Dominated Convergence Theorem, the following upper bound is used: 
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(85) 


26 


























where, as indicated in Veillette and Taqqu (2013), we have applied the inequality | exp(iz) — 

1 — z\ < for z € M. The right-hand side of (|85j) is continuous, for 0 < u < oo, and from 

2 

Lemma 4.1 of Veillette and Taqqu (2013) with G^^ d Jx) = Yk=\ 1 , keeping in 

mind the asymptotic order of eigenvalues of operator IC a (see, for example, Theorem 3.1 (i) 
by Leonenko, Ruiz-Medina and Taqqu, 2014), we obtain 

y2 

u G\(a/d) (exp(—u/2)) ~ uex.p(-u/2\ 1 (s$5)), as u —> oo 
uG X ( a /d) (exp(—u/2)) ~ [c(d,a)\D\ 1 - a/d } l/1 - a / d -^— 

xT (-—— ] (1 — exp(— u/2))~ 1 ^ 1 ~ a ^ d ' > ~ Cu~ 1 ~ a / d as u—> 0, (86) 

\ 1 — a/dJ 


for some constant C. Since 0 < Ya/d < 1, the right-hand side of (|86|h which does not 
depend on M, is integrable on (0,oo). Hence, by the Dominated Convergence Theorem, 


E 


ex 
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= exp 
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OO 

oo 
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exp(i#S£j) 


(exp(i6»u) - 1 - i Qu) —G x{a/d) (exp(-«/2)) 


du 


(87) 


which proves that equations (1811) and (I82[) hold. Equation (1831) follows, in a similar way 
to the proof of Theorem 5.1 (i) in Leonenko, Ruiz-Medina and Taqqu (2014), considering 
the expression obtained by the Levy density q in equation 


From the above equations, in a similar way as in Theorem 5.1 (ii)-(iv) by Leonenko, 

Y 2 

Ruiz-Medina and Taqqu (2014), it can be seen that S$S E TT>(M) is selfdecomposable. 

y 2 

Hence, it has a bounded density. It can also be showed that S$S is in the Thorin class 
with Thorin measure 

OO 

U{dx) = ^Y 5 _^( x )’ 

Z k= 1 2 \ k (S*£) 

x 2 

where 5 a (x) is the Dirac delta-function at point a. Finally, SSS admits the integral repre¬ 
sentation 

/ OO / \ p OO 

exp (-u) d I Y, ^k{S™)A > ' k ' , (u) J = / exp (—u)dZ(u), (88) 

\fc=i J 

where 

oo 

m = LA (S^)A^(t), t> 0, (89) 

k =1 

with A^ k \ k > 1, being independent copies of a Levy process. 
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